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Abstract 

In this paper, a thermodynamic analysis of Bouc- Wen models endowed with both strength and 
stiffness degradation is provided. It is based on the relationship between the flow rules of these 
models and those of the endochronic plasticity theory with damage, discussed in a companion paper 
(Erlicher and Point, 2007). Using the theoretical framework of that extended endochronic theory, it is 
shown that an elastic Bouc- Wen model with damage, i.e. without plastic strains, can be formulated. 
Moreover, a proper definition of the dissipated energy of these Bouc- Wen models with degradation 
is given and some thermodynamic constraints on the parameters defining the models behavior are 
emphasized and discussed. In particular, some properties of the energetic linear stiffness degradation 
rule as well as the so-called pivot rule, well-known in the seismic engineering field, are illustrated and 
commented upon. An improved energetic stiffness degradation rule and a new stiffness degradation 
rule are proposed. 
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1 Introduction 



Boud ( 119671 . Il97ll ) developed an univariate model for structural and ferro-magn etic applica- 



tions. The Bouc's formulation was similar to the one proposed by IVolterral (119281 ) to represent 



hereditary phenomena. However, the clock-time in the Volterra-Stieltjes integrals was substi- 
tuted by an internal or intrinsic ti me. 
several authors; see, among others 



Wen 



le Bouc model was modified by the contributions of 

1 I ' 



(1976) 



Karray and Boud (119891 ) and 



Casciatil (119891 ). 



leading to a more general class of models, hereafter named Bouc- Wen-type (BW) models. A 
further important modification concerned the introduction of the so-c alled strength and stiff - 



ness degradation effects, by means of suitable degradation functions (IBaber and Wen 



1981( 1 



In fact, collapse assessment in earthquake engi neering requires hy steretic models that include 



strength and stiffness deterioration properties ( jlbarra et al. 



20051 ). Some general assumptions 



on the form of these functions were provided, e.g. they have to be positive and monotone, but 



no other limitation was mentioned. The resultin g hysteresis models we re extensiv e 



seismic analysis of structural components, e.g. 



Reinhorn et al 



(119951 ): 



y use d in the 



Folientel (119951 ). More- 



over, the simplicity of BW models, due to the absence of an elastic domain, allowed their use in 
stochastic structural analysis. More recently, the BW models were used in applications of struc- 



tural contrc 


1, in particular in the r 


(Sain et al.. 


1997: 


Jansen and Dvke 



20001 ) or other kinds of damping devices (jlkhouane et al. 



20051 ). Moreover, several contributions of the last three decades were concerned w ith the iden- 



tificat ion of BW model parameters, for applications in seismic engineering; e.g. 



J200J). 



Masri et al. 



Valanid (119711 ) proposed 



Starting from a different viewpoint and independently from Bouc, 
the endochronic (EC) theory of visco-plasticity, which postulates the existence of an intrinsic 
time governing the rate-independent evolution of stress and strain in materials, whereas the 



Newtonian time is used to model the viscous behavior. In the case of plasticity without viscous 



3 



effects, it was proven that flow rules can be formulated by using th e mathematical formalism 



of pseudo-potentials. This was done for the standard EC theory by 



Erlicher and Pointl (12 006) 



and for the EC theory with isotropic damage in the companion paper (jErlicher and Point 



20071 ) 



A major question concerning the BW models is their thermodynamic admissibility. Sev- 
eral criticisms were moved to these models as their formulation was not based on a p roper 



thermodynamic approach. Re cently, this problem was analyzed by 



Ahmadi et al. 



(119971 ) and 



Capecchi and De Felice! ( 120011 ) for models without degradation terms. On the other hand, the 



typical constitutive laws of EC models are characterized by the absence of an elastic domain 



and the corresponding hysteresis loops are smoo t 



reason several authors, among others, 



i and open , like in the BW models. F or thi s 



Casciatil (119891 ) and 



Sivaselvan and Reinhornl (120001 ). 



observed that a relationship must exist between EC and BW models. A formal proof that all 
the Bouc-Wen type models with a strength degradation term admit an equivalent endochronic 



formulation was provided by 



Erlicher and Point 



(120041 ) . This was sufficient to prove that BW 



models with strength degradation are thermodynamically admissible, in the sense that they 
fulfill the second principle. 

Nonetheless, a problem remains unsolved, viz. the formal proof of the thermodynamic admis- 
sibility of BW models endowed with a stiffness degradation term: this is th e subject of the 



presen t paper. In this respect, the endochronic theory with damage discussed in 



Erlicher and Point 



( 120071 ) is used: it is proved that for every BW model with strength and stiffness degradation, 
an equivalent EC model with damage can be defined, i.e. a model exhibiting the same flow 
rules. Therefore, the BW model is thermodynamically admissible if the corresponding EC 
model fulfills the second principle. Exploiting the aforementioned equivalence, the following 
results are derived too: (i) a rigorous thermodynamic definition of the plastic strain and of 
the dissipated energy for BW models exhibiting strength and stiffness degradation; (ii) a new 
elasto-damaged BW model; (iii) a bound for the stiffness degradation rate of BW models, 
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proving that the traditional semi-empirical stiffness degradation laws, like the energetic lin- 
ear rule or the so-called pivot rule, may lead to some pathological thermodynamic behavior 
(these situations are illustrated by several numerical examples); and (iv) an improved energetic 
stiffness degradation rule and a new stiffness degradation rule. 

After the introduction in the first Section, the existing BW models are reviewed in the first 
three parts of the second Section, while a new class of stress-strain BW models with stiffness 
degradation is defined in the fourth part. The equivalence between the flow rules of BW models 
with both degradations and the EC theory with damage is discussed in the third Section. In 
the fourth Section, single degree-of-freedom (1-DoF) and two-degrees-of-freedom (2-DoF) BW 
structural models are analyzed in detail, using the theoretical tools previously developed. This 
analysis is supplemented by some numerical examples and by an application to partial-strength 
beam-to-column steel joints. 



2 Bouc- Wen-type models 

A classification of Bouc- Wen- type models is suggested in Figure 1. The main distinction con- 
cerns non-degrading models (ND-BW models) and degrading ones. Moreover, among degrading 
models, the further differentiation between D-BW models having only strength degradation 
and models characterized by strength and stiffness degradation or by stiffness degradation only 
(DD-BW models) is pointed out. Then, the models formulated for describing the material be- 
havior by stress-strain laws and the structural behavior by generalized force-displacement 
laws are distinguished. Note that the box corresponding to stress-strain DD-BW models rep- 
resents a new class of models, defined by equation (20). In this paper, the equi v alenc e between 
EC models with isotropic damage (DD-EC), as defined in lErlicher and Pointl (120071 ). and the 
stress-strain DD-BW models is proven and then used to provide a thermodynamically well- 
posed formulation for the DD-BW models. 
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The stress-strain models described in this paper are presented only for a plastically incom- 
pressible behavior, since it is the most important case in the literature concerning stress-strain 
BW models. The treatment of the more general case, with a non-zero hydrostatic plastic flow, 
is conceptually similar, but requires a more complex formalism not needed here. 



2. 1 ND-BW models for structures 



Among the different models of hysteresis proposed by iBoud (119711 ). the simplest one can be 
formulated by a Stieltjes integral as follows: 



w (t) = A u(t)+Z (t) 



(1) 



where A is a non- negative constant; u and w are two time-dependent functions, which are 
considered as input and output functions, respectively. In structural engineering applications 
the input usually has the meaning of a generalized displacement, while the output w plays the 
role of a generalized force, defined as the sum of a linear term A u(t) and a hysteretic term 
Z(t). The integral in (1)2 depends on the time-function £ s (t), which is named internal time 
and is assumed to be positive and non- decreasing. The function /x, called the hereditary kernel, 
takes into account hysteretic phenomena. One of the definitions of £ s proposed by Bouc is the 
total variation of u: 





du 


f 


dtf 



dt' or, equivalently, £ s 



with £ s (0) = 0. 



(2) 



where the superposed dot indicates time-differentiation. (2) i mplies the ra te-independence 



of £ s and as a result, Z and w are in turn rate-independent. 



Boud (|197lD defined /x as 81 



continuous, bounded, positive and non- increasing function on the interval £ s > 0, having a 
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bounded integral. In the special case of an exponential kernel 



with A,p s >0 



(3) 



a differential formulation of (1) can be easily deduced. Then, (2) and (3) imply: 



w = A u + Z 

Z = Au-(3 s Zi 



(4) 



with £ s = \it\ 



One can observe that for an initial value in the interva 
teretic force Z remains in the same interval, 
of (4): 

w = A u + Z 



[-Z u , Z u ), with Z u = A/P s , the hys- 



Boud (119671 ) proposed a more general formulation 



while 



(5) 



Z = A it — (3 a Z — j s \Z\ u 



with 7 S < (3 S 



Wen! ( 119761 ) suggested a further modification by introducing a positive exponent n: 



w = A u+ Z 

Z = A u — (P s sgn [Z it) + j a ) \Z\ n u 



(6) 



where sgn(-) is the signum function. Wen did not impose any condition on the value of 7 S and 
assumed integer values for n; however, all real positive values of n are admissible. When n is 
large enough, force-displacement curves similar to those of an elastic-perfectly-plastic model 
with the additional linear term Aqu are obtained. Provided that /3 S + 7 S > 0, the limit strength 
value Z u of (6) becomes: 

Z n =(-^—Y (7) 



Ps + ls 



The parameter p a is positi ve by assumption, w 



from the condition £ s > ( jErlicher and Point 



rile t he admissible values for ^ s can be derived 



20041 1 . See also Table 1, in this respect. 
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2.2 Plastically incompressible ND-BW stress-strain models 



In order to link t he cla ssical plasticity theory and the Bouc-Wen model described in (6), 



Karray and Boud (119891 ) proposed a tensorial generalization of (6) for isotropic and plastically 



incompressible materials: 



tr (cr) = 3E t tr (e) + 3K tr (e) 



dev (cr) = 2E 2 dev (e) + z 



2G dev (i)-(3z£ 



with (3 > and the following special choice for the intrinsic time 
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(8) 



£ = ( 1 + — sgn (z :dev (e)) J |z : dev (i) 



in-2 



(9) 



The terms proportional to i?i > and E 2 > play the same role as A . £ and cr are the 
second-order symmetric strain and stress tensors, respectively; tr and dev are the trace and 
the deviatoric operator, respectively; z is a traceless tensor defining the hysteretic part of the 
stress, while ||-|| is the standard L 2 — norm of a symmetric second-order tensor. K and G are 
the bulk modulus and the shear modulus, respectively. Note that the stress-str ain relationship 



for the traces is linear, according to the assumption of plastic incompressibility. ICasciatil (119891 ) 



proposed a stochastic analysis of the same model. For this reason, we name it as Karray-Bouc- 
Casciati (KBC) model. The norm of the tensor z (t) is bounded as follows: 

' 2G \" 



\dev (cr) — 2E 2 dev (e)|| < a u 



(10) 



for t > 0, provided that ||z (0)|| < a u . This inequality shows that a limit strength value exists 
for the KBC model and only regards the deviatoric part of the stress cr, consistently with the 
plastic incompressibility requirement. A thermodynamic formulation of this model, based on 
a suited definition of the Helmholtz free energy and of a pseudo-potential, was proposed in 
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Erlicher and Pointl (120061 ). 



2.3 DD-BW models for structures 



The BW model (6) was modified by lBaber and Wen! (Il98ll ). who introduced the positive and 
increasing functions r] s and v s , both having a unit initial value: 



w = A u + Z 

Z = ±-[Aii- v s ((3 s sgn (Z u) + j B ) \Z\ n u] 



(11) 



where v s represents a strength degradation effect, while r\ s is associated with a stiffness degra- 
dation effect. According to Figure 1, an attentive reader can see that this model belongs to 
the class of DD-BW models for structures. A degradation function similar to rj s could also be 
associated with the post-yielding stiffness Aq. However, in order t o simplify the present a nal- 



Bab er and Wen 



ysis, A is supposed to be constant. In the original formulation of 
and v s are defined as linearly increasing functions of e s , the energy dissipated by the hysteretic 
model: 

v s = 1 + c UtS e s , rj s = 1 + c rus e s (12) 

with Cj, jS > and c V)S > 0. The authors did not provide any definition of e s on a thermodynamic 
basis. This topic will be discussed later, in the fourth Section. The case of increasing strength, 
viz. < v s < 1, is also admissible and it is analogous to isotropic hardening in classical 
plasticity models. Provided that (3 S + 7 S > 0, the maximum strength value, modified by the 
strength coefficient r, becomes 



.4 



J u,red 



r Z„ 



(13) 



Several applicat i ons o 



others, iFolientd (1l995l ) and 



s i> s {j3 s + 7 S ) / 

' this version of the BW model can be found in the literature; see, among 



Sivaselvan and Reinhornl (120001 ) . The latters proposed a formula- 
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tion of (11) entailing a clear physical meaning of the parameters. Moreover, the following 
degradation rules were suggested: 



§2 e 3 
-P2 e ui 



(14) 



rj s = rj s (u, w) 



R k (A+A )-A 



m + a Z u,red 
{A+A )u + a Z u red 



(A+A )-A 



where u max = rnax < r <i (\u (r)|) is the actual value of the maximum displacement modulus; 
u u it and e u u are defined as the ultimate displacement and the ultimate dissipated energy, 
corresponding to failure; (3i and /3 2 are positive parameters related to strength degradation; 
a is a positive parameter related to the stiffness degradation rule (14)2, called the pivot rule. 



(A+A ) 



is related to the position of the actual force- 



The sign of Z u , red = Z u , red sgn (u 
displacement point (u, w) with respect to the line representing the initial elastic behavior 
w = (A + A )u. 

In summary, the general univariate (also called uniaxial or 1-DoF) BW model reads: 



w = A u + Z 
Z = ±(Au-f3 s Z u s i 

where 7] s , v s and £ s can vary from one model to another, as shown in Table 1. 



(15) 



n r 



A 2-DoF generalization of (15) was defined by iPark et al.l (119861 ) to represent the behavior of 
a system constituted of a single mass m subjected to an excitation acting in two orthogonal di- 
rections. A simple model associated with this situation is illustrated in Figure 2a. For instance, 



this model is suited to reproduce t he geometri c ally 



loaded reinforced concrete column ( 



laminated rubber bearing (lAbe et al 



Park et al. 



inear uncoupled behavior of a bi-axially 



19861 ) or the uncoupled biaxial behavior of a 



2004]). In the simplest case where the stiffness and the 



10 



strength are the same in both directions, i.e. isotropic behavior, one has 



W = A V + Z 



Z=±(AU-p a Zu.t) 



(16) 



where W = (w x ,w y ) is the total force vector, U = (u x ,u y ) is the structural displacement and 
Z = (Z x , Zy) is the hysteretic force; the parameters A$ and A have the dimension of a stiffness; 
77 s , v s are scalars, since the damage and strength evolution are assumed to be equal in both 
directions. The intrinsic time £ s = £ s (Z x , ii x , Z y , u y ) and the degradation functions u s , rj s may 
be chosen according to one of the definitions given in Table 2. 

An equivalent way of writing (16) is 

w x = A u x + Z x 



Wy = A Uy + Zy 

Z x = ±- s \A u x - f3 s v s Z x t 

Zy = ^ (yA Uy — (3 S U S Zy^ S 



(17) 



A special case of (17) is given by (IPark et al 

w x = A u x + Z x 

Wy = A Q Uy + Zy 



19861 ) 



;is) 



Z x = ± {A u x - v s [f3 s Z x \Z x u x \ + 7 S Zlu x + (3 S Z x \Z y u y \ + 7 S Z x Z y u y \) 



Zy= ±(A.Uy— V s (3 S Zy \ZyUy\ + 7s Z^Uy + /3 S ^ | Z X U X \ + 7s ^X^: 



which correspo nds to a pa r ticula r choice of the intrinsic time £ s , as indicated in the second 



row of Table 2. IPark et al. 



1986T 1 assumed f3 s = 7 S with A and 7 S having a special coupled 
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evolution, introduced to obtain force-displacement curves confined in a given envelope. How- 
ever, we shall assume that A and 7 S are constant. The strength degradation u s was defined 
as in (12)i, while the stiffness degradation i] s was given by: 



Vs = Vs (Vax,x, WmaxJ = J Vl + Vy = \ + > 1 (19) 

V \ \ Uyield ) \ Uyi e i d ) 

where u y i e u is the yielding displacement, u mSuX , x and u maX! y are the actual maximum values of 
the displacement in the two directions x and y. The stiffness degradation is therefore governed 
by the "uniaxial" ductility ratios r\ x and r] y , while i] s is a "bi-axial" ductility ratio. 



2.4 Plastically incompressible DD-BW stress-strain models 



The KBC model (8)-(9) is a stress-strain BW model characterized by the absence of strength 
and stiffness degradation. In this subsection, a direct generalization of (8) is proposed, leading 
a DD-BW model: 

tr (a) = 3E 1 tr (s) + ±3K tr (e) 

( dev (a) = 2E 2 dev (e) + z ( 20 ) 

z = i (2G dev {i)-(3zvi) 
The linearity of the hydrostatic stress-strain law is due to the assumption of plastic incom- 
pressibility. The thermodynamic admissibility of the DD-BW model defined in (20), i.e. the 
conditions under which it fulfills the second principle of thermodynamics, is discussed in the 
two following Sections. Some simple degradation functions may be defined by analogy with 
the structural energetic linear rules discussed in the previous subsection, viz. 

v = 1 + c u e , r] = 1 + c v e (21) 

where e is the dissipated energy per unit volume and the coefficients c u , c v are analogous to 
c^s, c VjS introduced in (12). The expressions of rj, v and £ for different models are collected in 
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Table 3. As already stated for structural models, the case of increasing strength, viz. < v < 1, 
is also admissible and it is analogous to isotropic hardening in classical plasticity. 

Models with strength degradation lead to stress-strain laws with softening and, if applied in 
the continuum setting, to strain localization. It is well-known that an objective description of 
softening materials and the resulting localized failure modes require special attention. Thus, 
some of the model parameters must be considered not as pure material properties but as 
discretization-dependent or, alternatively, the underlying theory must be enriched by special 
terms, like weighted spatia l averages or higher-order gradients, acting as localization limiters 



( iJirasek and Bazantl . 120021 . Chap. 21). These aspects are not the main subject of this paper 



and therefore are not treated here. 



3 Equivalence between DD-BW stress-strain models and endochronic models 
with damage 



he endochronic models w ith isotropic damage (DD-EC) are discussed in a companion paper 



( lErlicher and Point 



20071 ) . We are interested here in the analysis of a DD-EC stress-strain 
law having the same format as the force-displacement rules (15) and (16), which is defined by 
the following relationships 



a = E : e + (1 - D) C : (e - e p ) 



(22) 



tr (i p 



and 



1 devja-E-.e) ( 
1-D 2G//3 g 



The non-degrading linear term E : £ is analogous to the term Aqu characterizing the models 
(15) and (16). The tensor £ p is the plastic strain. The trace of the plastic strain flow is zero 
because of the plastic incompressibility assumption; C = (k — 1 <g> 1+2GI is the standard 
elasticity fourth-order tensor for isotropic materials; E = {Ei — ^l 2 -) 1 ® 1 + 2E 2 1 is a fourth- 
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order tensor defining an additional elastic stiffness; 1 is the second-order identity tensor; I is 
the fourth order identity tensor and £g> represents the tensor product; D is the scalar damage 
variable, supposed to be non-negative and less than one; ( is called the i ntrinsic time measure 



and C > has the role of plastic multi plier ( 



Er 



called the hardening- softening junction (IBazant 



icher and Point 



20071 ); g > is sometimes 



19781 ). (22) is equivalent to 



tr O) = 3^ tr (e) + (1 - D) 3K tr (e) 



dev (<x) = 2E 2 dev (e) + z 



1 — D) 2G dev (e) - (3 z £ - D ^ 



(23) 



where the time derivative of the hysteretic deviatoric tensor z is explicitly written. The com- 
parison of (20) and (23), in particular the definition of z, shows that the two models, DD-BW 
and DD-EC, are equivalent if 

1 : C 



1 

V — — 

9 



V 



D 



9 



(24) 



1-D 1 * 1-D P{l-Df 

with v (0) = 1, D (0) = and £ (0) = £ (0) = 0. (24) relate the functions D, g and ( associated 
with the DD-EC model and the functions rj, v and £ associated with the DD-BW model. Obvi- 
ously, this equivalence assumes that the evolution of D, g and ( can be defined in a way that is 
general enough to represent all possible evolution of the given 7], v and £. To be more precise, if 
for instance t] and v are proportional to the dissipated energy, which for DD-BW models may 
increase during unloading phases, then it should be possible to define g and D of the equiva- 
lent endochronic model such that they depend on the dissipated energy, with possible non-zero 
increments also during unloading phases. The proper thermodynamic fra mework for DD-EC 



mode ls allowing this general behavior is presented in a companion paper (jErlicher and Point 



20071 ) . (24) i shows that the function v defining the strength increment (degradation) for BW 
models is strictly related to the function g defining the isotropic hardening (softening) for en- 
dochronic models. (24)2 shows a similar equivalence between the stiffness degradation function 
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of DD-BW models and the isotropic damage variable for DD-EC models. The interpretation 
of (24)3 is more involved. It states that the intrinsic time flow £ of a DD-BW model can be 
expressed as the sum of: (i) a contribution related to the plastic multiplier ( of the equivalent 
DD-EC model; (ii) a part related to damage flow D of the equivalent DD-EC model. Let us 
set (24)3 in the alternative form 

C = - - IT" (25) 
77 pr\v 

One can see that (25) is the definition of the plastic multiplier ( of the DD-EC model equivalent 
to a given DD-BW model, as a function of (3, £, rj and v. Briefly, (25) defines the plastic 
multiplier of stress- strain DD-BW models. When ( = 0, there is no plastic flow, viz. an elastic 
with damage behavior is retrieved. Moreover, ( = if and only if 

V = 9 = 0vi (26) 

(1 — D) y 1 

i.e. for DD-BW models, it is possible to obtain zero plastic strains, provided that the special 
stiffness degradation rule (26) is adopted. Figure 3 shows this behavior, with the parameters 
El = E2 = 0, 2G = 24000 MPa, K = 20000 MPa, £ is defined by (9) with n = 1.2, 



(3 = 81.2409 MPa 1 - 11 , 7//? = 0.5 (entailing a u = 100^(2/3) MPa), v = l + c u e with c„ = 0.05 
m 3 MJ~ 1 , e is given in (33) underafter and 77 is defined by (26). 



4 Thermodynamic analysis of DD-BW models 



4-1 DD-BW stress-strain models 



The ther modynamic analysis of the endochronic model with isotropic damage (DD-EC) pre- 



sented by 



Erlicher and Point 



(120071 ) led to the following conditions ensuring thermodynamic 
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admissibility: 

D > : the damage is non- decreasing 

(27) 

£ > : the plastic multiplier is non-negative 
These conditions derive from the assumption that the dissipated energy increments owing 
to plasticity and to damage are separately non-negative. It is well known that this assump- 
tion is not needed, since only the total dissipated energy must be non- decreasing. However, 
this assumption suffices to defi ne models fulfilling the second principle; it is usually adopted 



( jLemaitre and Chaboche 



19901 ) and is also used hereafter. (24) 2 and (27) i lead to 

f] > : the stiffness degradation function is non- decreasing (28) 

while (24) 3 and (27)i_ 2 entail 

£ > : the intrinsic time of BW models is non-decreasing (29) 

For DD-BW models, an expression of £ fulfilling (29) is usually given, while ( is unknown. It 
follows that the positivity of (, recall (27)2, must be ensured by a relevant condition on the 
stiffness degradation rule, which can easily be derived using (25): 

V < v i : upper bound for the stiffness degradation rate (^0) 

In summary, two conditions on the stiffness degradation rule are provided: the first one, i.e. 
(28), is trivial, while the second one, (30), is less obvious. They are obtained using (24) and 
(25), that relate a DD-BW model, defined by £, r] and is, with the equivalent DD-EC model, 
defined by (, D and g. 

The inequality (30) suggests the definition of a new stiffness degradation rule (recall that 
T]>1): 

V = cJ^) Pu£ (31) 
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with cp E [0, 1] and m > 0. For m = 0, (31) entails (see also (24)-(25)) 

r 1 = cpf3vi ^ (=(1-D) ai-cp) (32) 

i.e. the plastic multiplier is proportional to (1 — D): with the increase of damage, the plastic 
multiplier decreases. In other words, this rule postulates that plastic strains are larger when 
the material is slightly damaged and viceversa. When also cp = 1 holds, then ( = and an 
elastic with damage behavior is retrieved, as shown in Figure 3. 

Moreover, one can observe that, still using (24)-(25) and recalling that z =dev (cr — E : e), the 
energies dissipated by plasticity and by damage read: 



e = e p + e D 



with 



(33) 



e D = Y d D = \{tr - E : e) : C 



-l . 



■in- F • p\ j> _ If*:*! Ma-Ee)) 2 \ ■ 
■ [(T ~ ■ £ > (T^Df _ 2 [2G + 9K ) V 



Erlicher and Point 



(120071 ) . Y d is 



See also the corresponding definitions for DD-EC models in 
the dissipative part of the thermodynamic force dual of the damage variable D. Note that 
the increments of the energy dissipated by plasticity e p depend not only on the BW intrinsic 
time flow £ but also on the flow f] of the stiffness degradation function. Figure 4 shows the 
evolution of the two terms of the dissipated energy as a function of time. The parameter values 
are the same of those of Figure 3, except for the ratio j//3 — —0.5 and j3 = 243.7227 MPa l ~ n \ 
v = 1 + c u e with c v = 0.05 m 3 MJ _1 ; 77 = 1 + c v e p with c n = 0.05 m 3 MJ _1 . Note that there 
are non-zero plastic strain increments during unloading pha ses. Moreover, damage inc rements 



are always non-zero, also during unloading phases; see also 
respect. 



Erlicher and Point 



(120071 ) in this 
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4-2 DD-BW models for structures 



Comparing (20) (that is equivalent to (22) under the conditions previously discussed) with 
(15) and (16), one can see that structural models have the same format as stress-strain models, 
provided that stress is replaced by force and strain by displacement, respectively. Analogous 
substitutions are made for the stiffness tensor, the (3, 7 parameters, the degradation functions 



and the intrinsic time (see the first two co 
device subjected to bi-directional shear, see 



umns of 



Abeetal 



' able 4 ). In the ideal case of a bearing 



(120041 ) and Figure 2b, the relationships 



between the stress-strain and force-displacement parameters are given in the third column of 
Table 4. In general, if a direct derivation of structural models from stress-strain models is 
not possible, the formal analogy between stress-strain and force-displacement rules must be 
postulated. Accounting for this analogy and referring to the isotropic 2-DoF model (16), we 
assume a force- displacement relationship similar to the stress-strain law (22), that reads: 



W : = A XJ + Z = A U-\ — A (U - IP) . 

Vs 



(34) 



The seventh line of Table 4, with the definition of D s , has also been used. Moreover, from 
(16), one has 



W = A XJ + — A\J - (3 S Z v s £ 
Vs 1 



(35) 



Equation (35) and the time-derivative of (34) lead to the definition of the plastic displacement 
flow 



Z Z 

U p = — (flsVsis - Vs) = -rflsVsVs — ~ 

A v ' A \t] s f3 s u s r] s 



A 



PsVsVsCs 



where ( s > is a structural plastic multiplier; compare with (25). By analogy with (33) and 
using Table 4, the rate of energy dissipated by the force-displacement model can be defined 
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as follows: 



&s,p ~r &s,D 



with 



e 



s,p — 




Z T -Z 

A 



(A 



's 



f s6 - Vs 



) 



(36) 



z T -z £) 3 

2A (1-D S ) 




All the previous quantities are defined for a 2-DoF BW model; 1-DoF models correspond to 
deal with scalar quantities instead of vectors. Concerning 2-DoF models, it is interesting to 
analyze the relationship between the case of an intrinsic time of the Park- Ang- Wen (PAW) 
type (see the second row of Table 2), versus a 2-DoF model having an intrinsic time of the 
KBC type with n = 2 (see the last row of Table 2). It is easy to prove that the intrinsic 
time flows £ s of the two models become identical in the case of proportional loading. However, 
they are different in the case of non-proportional loading, leading to different hysteretic loops, 
as illustrated in Figure 5. The parameters chosen for the numerical simulation are n = 2, 
A = 0, A = 3.571 kN/mm, (3 S = 0.0119 kN^mm' 1 , -f s = -/3 s /4. Moreover, the strength 
degradation is defined by (12) x with c UyS = 0.0001 J -1 , while the stiffness degradation by (19) 
with Uyi e id = 11.2 mm. 

4-2.1 Discussion about the stiffness degradation rule 

Conditions (28) and (30) hold at the stress-strain law level. However, taking into account the 
above-mentioned analysis, the same constraints are postulated at the force-displacement level. 
As a result, one has 



Likewise, also the stiffness degradation rule (31) can be extended to BW models for structures: 



According to the proposed analysis, every stiffness degradation rule defined for structural BW 



Vs > , f] s < (3 S u s £ s 



(37) 




with C/3 G [0, 1] and m > 



(38) 
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models should be checked with respect to (37). As a result, some bounds for the parameters 
characterizing these degradation rules may be found. Two examples are considered here: the 
energetic linear rule, where a possible violation of (37)2 is highlighted, and the pivot rule, 
where (37) i may not be fulfilled. 



Baber and Wen 



(Il98ll ). postulates 



The energetic linear rule defined in (12)2, first adopted in 
that the stiffness degradation function rj s depends linearly on the dissipated energy e s 
e s ,p + e s ,D, defined in (36): 



Cti, s e s Cr?,s (&s,p &s,d) ¥ ^js ^ Z T -Z P s ^ s ^ s (39) 

1 + C^ jS 2~ — 

As a result, (37)2 imposes that the coefficient multiplying (3 s v s ^ s in (39)2 be less or equal than 
one, entailing c^s Z T -Z < 2 A. This is an upper bound for the parameter c^s > 0, varying with 
the actual value of Z. When it is violated, (37)2 is not satisfied, leading to a negative plastic 
multiplier, i.e. ( s < 0, and to negative increments of the energy dissipated by plasticity, i.e. 
e S)P < 0. This situation is considered as non-admissible, since it is assumed that the energies 
dissipated by plasticity and damage must be separately increasing. A numerical simulation 
where this pathological situation occurs is presented in Figures 6a-c, where a 2-DoF BW 
model with intrinsic time of KBC type is considered; in detail, A = 0, A = 3.5714 kN/mm, 
f3 s = 0.1308 kN 1 ~ n mm~ 1 , 7 S = — /3 s /4 and n = 1.2. The strength degradation is governed by 
(12)x with c UtS = 0.005 J -1 , while the stiffness degradation is given by (12) 2 with c V:S = 0.1 J" 1 . 
Conversely, defining a stiffness degradation rule that depends linearly on the energy dissipated 
by plasticity 



' ItCm.s a — 



(40) 



the limit condition (37)2 is always fulfilled. Hence, for a stiffness degradation proportional 
to the plastic dissipated energy, no limitations are needed on c V)S . The numerical example 
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illustrated in Figures 6d-f confirms this result. All the model parameters are the same used 
for the previous case. 

The pivot rule (14)? is sometimes used for modeling the cyclic behavior of structures; see 



e.g. 



Sivaselvan and Reinhornl (120001 ) in this respect. Some problems associated with this rule, 



occurring when it is applied with displa c emen ts histories with cycles of decreasing amplitude, 



were pointed out by lWang and Folientd (120011 ). A BW model is considered here, with A = 0, 



A = 50, 7s = P s = 0.0189, n = 1.2 (entailing Z u = 400) and v s = 1 + c U)S e s with c u ^ s = le — 6. 
Th e parameter a of the pivot rule is equal to 5. In this example, the same difficulties mentioned 



by IWang and Folientd (j200ll ) appear, viz. the stiffness increases when the cycle amplitude 
decreases: see Figure 7c. Moreover, a deeper problem is also highlighted: even if the energy 
dissipated by plasticity and the total dissipated energy are non-decreasing, and the second 
principle is fulfilled, see Figure 7d, this rule may entail negative damage increments, both 
during increasing and decreasing amplitude loops. As a result, D s < r) s < with a 
violation of (37) i, as shown in Figures 7e,f. 

Let us now compare the pivot rule to the stiffness degradation rule (38). Figure 7b shows that it 
is possible to choose the parameters for both rules in such a way to have almost identical force- 
displacement loops during an initial phase characterized by increasing displacement amplitude 
cycles. In the specific case considered in the Figure, the new stiffness degradation rule (38) 
is characterized by eg = 0.06 and m = 1.3. The damage evolution during this phase for 
both degradation laws is depicted in Figure 7e: the new degradation rule entails a monotonic 
increase of damage, consistent with the assumption that dissipated energies by plasticity and 
damage separately increase; while the pivot rule is characterized by phases with alternate 
increase and decrease of damage. Then, decreasing amplitude loops are considered for both 
models and they provide very different behaviors, as illustrated in Figure 7f: the rule (38) still 
entails an increase of damage, while the pivot rule induces an averaged decrease of damage, 
corresponding to an averaged increase of stiffness. A similar comparison can be made between 
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the pivot rule and (40). Figure 8 shows the numerical results with c VjS = 8.4e — 6. The same 
remarks made for Figure 7 hold: Figure 8d shows the evolution of the total dissipated energy 
for the case of the pivot rule (thin line), the rule (40) (thick line) and the rule (38) considered 
in the previous example (dotted line). As expected, the pivot rule has a different evolution in 
the last part of the displacement history, where the stiffness increases instead of decreasing 
like for the other two rules. 

An application application of the BW model (11) with the stiffness degradation (38) and the 
strength degradation v s — ( 1 — y-f^) Ps j , inspired to the rule (14) i is shown in Figure 9. In 
det ail, the exp e rimen tal loops refer to a partial-strength beam-to-column steel joint, studied 



by 



Bursietal 



(120021 ). The assumed parameter values are A = 1.0 kNm/mrad, A = 21.0 
kNm/mrad, (3 S = 0.4888 (kNm) 1 ~ n (mrad)~\ 7s = 0, n = 0.8, Z u = 110 kNm, e uH = 210000 
Nm, (3% = 0.35, op = 0.02 and m = 1.8. For the pivot rule, we assumed a = 4. An attentive 
reader can observe that a good agreement between experimental (Figure 9b) and numerical 
(Figure 9c, d) hysteretic loops is obtained. In addition, the evolution of damage predicted by 
the model with the rule (38) and with the pivot rule is illustrated in Figure 9e. Again, the 
drawback of the pivot rule is evident. In summary, one can state that (38) and (40) may 
be used as a valid alternative to the pivot rule to predict accurate loops; furthermore, they 
represent a significant improvement in terms of thermodynamic admissibility. 



5 Conclusions 



An exhaustive classification of Bouc- Wen-type models was provided in this paper. The mod- 
els endowed with both strength and stiffness degradation were compared with an extended 
endochronic theory with isotropic damage. The thermodyn amic formulation of t his ex tended 



endochronic theory was discussed in a companion paper (jErlicher and Point 



20071 ). using 



pseudo-potentials depending on state variables and on parameters related to the past history 
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of the material. Using this theoretical framework, a rigorous thermodynamic definition of the 
energy dissipated by plasticity and damage was provided for Bouc-Wen models having both 
strength and stiffness degradation, as well as the definition of an elastic with damage Bouc- 
Wen model. Moreover, an important constraint on the stiffness degradation rules for Bouc-Wen 
models was highlighted and its consequences on the energetic linear rule were investigated. 
The non-monotonic damage evolution characterizing the pivot rule was illustrated. Finally, 
an improved formulation of the energetic linear rule and a new stiffness degradation rule were 
proposed; their effectiveness with respect to the pivot rule was proved by numerical examples 
and an application to beam-to-column steel joints. 

Acknowledgements 

The authors would like to thank prof. Nelly Point and acknowledge the several discussions 
which helped the improvement of this work. The authors are also grateful to a reviewer for 
his constructive criticism. 

6 Appendix: Notations 

The following symbols are used in this paper: 
Aq = additional stiffness for structural BW models 
A = initial stiffness for structural BW models 
C = fourth-order elasticity tensor 
D = internal variable associated with isotropic damage 
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E = additional fourth-order elasticity tensor 

Ei = bulk modulus of the additional elastic behavior 

E 2 = shear modulus of the additional elastic behavior 

K = bulk modulus 

G = shear modulus 

I = fourth-order identity tensor 

Y d = dissipative part of the thermodynamic force dual of the damage variable 
1 = second-order identity tensor 

e = total dissipated energy per unit volume 

e s = total structural dissipated energy 

e p = energy per unit volume dissipated by plasticity 

= structural energy dissipated by plasticity 
er> — energy per unit volume dissipated by damage 
z s ,d = structural energy dissipated by damage 
g = hardening-softening function for EC models 
u = total displacement for 1-DoF BW models 
U = total displacement vector for 2-DoF BW models 
w = total force for 1-DoF BW models 
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W = total force vector for 2-DoF BW models 

Z = hysteretic force for 1-DoF BW models 

z = hysteretic part of the stress tensor 

Z = hysteretic force vector for 2-DOF BW models 

a = coefficient associated with the pivot degradation rule 

(3 = coefficient characterizing the flow rules of BW or EC stress-strain models 

(3 S = coefficient characterizing the flow rules of structural BW models 

7 = additional coefficient defining the intrinsic time of BW or EC stress-strain models 

7 S = additional coefficient defining the intrinsic time of structural BW models 

e = total small strain tensor 

s p = plastic strain tensor 

(j, = hereditary kernel 

v = strength degradation function for BW stress-strain models 
v s = strength degradation function for BW structural models 
rj = stiffness degradation function for BW stress-strain models 
7] s = stiffness degradation function for BW structural models 
a = Cauchy stress tensor 
d = intrinsic time scale for EC models 
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£ = intrinsic time for BW stress-strain models 

£ s = intrinsic time for 1-DoF and 2-DoF BW models 

( = intrinsic time measure for EC models; plastic multiplier for stress-strain BW models 
( s = structural plastic multiplier for 1-DoF and 2-DoF BW models 
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List of Tables 



Summary of the 1-DoF BW models for structures. £ s > 0, provided that 
|7s| < As- For a thermodynamically admissible model, the stiffness degradation 
rule fulfills the condition fj s < f3 s p s ^ s (see equation (37)2). 

Summary of the 2-DoF BW models for structures. Both the conditions on the 
parameter 7 S and on the stiffness degradation rule recalled in the caption of 
Table 1 hold for 2-DoF models too. The last definition of £ s is named here of 
the KBC-type. 

Summary of the BW stress-strain models. Observe that the first BW model 



Valanisl fll97lh . £ > 0, 



coincides with the basic endochronic model proposed by 
provided that I7I < (3. For a thermodynamically admissible model, the stiffness 
degradation rule fulfills the condition 1) < (see equation (30)). 

Stress-strain vs. force-displacement rules. In the ideal case of a bearing device 
subjected to bi-directional shear, see Figure 2b, stresses and strains are 
uniform in the device volume. Both £ s and £ are postulated to be of the KBC 
type. i] s and v s are expressed as functions of global quantities like forces 
and displacements, while £, r\ and v depend on local quantities. 
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Table 1 

Summary of the 1-DoF BW models for structures. £ s > 0, provided that |7 S | < [3 S . For a thermo- 
dynamically admissible model, the stiffness degradation rule fulfills the condition f] s < /3 S ^ S ^ S (see 
equation (37)a). 
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Table 2 

Summary of the 2-DoF BW models for structures. Both the conditions on the parameter 7 S and on 
the stiffness degradation rule recalled in the caption of Table 1 hold for 2-DoF models too. The last 
definition of £ s is named here of the KBC-type. 
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Table 3 

Summary of the BW stress-strai n mode" 
endochronic model proposed by 



\\dev(i)\\ 



1 + 77 sgn (z :dev (e)) ) |z : dev (e)| ||z|| n 2 , n > 



1 + %. sgn (z :dev (e)) ) |z : dev (e)| ||z|| re 2 , n > 



1 + a s^n (z :deu (e)) ) |z : det> (e)| ||z|| n 2 , n > 



s. Observe that the first BW model coincides with the basic 



Valanid (|197ll ). ^ > 0, provided that I7I < f3. For a thermodynam- 
ically admissible model, the stiffness degradation rule fulfills the condition f\ < f3u^ (see equation 
(30)). 
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Stress-strain law (20) 


Force-displ. laws (15)-(16) 


Bearing shear device 


a- [Pa] 


W X ,Wy [N] 


W X = d XZ S, Zy = (Jy Z S 


z [Pa] 


Z X ,Zy [N] 


Z X = Z XZ S, Zy = Zy Z S 


e 


u x ,u y [m] 


U x — 2/iiS x z - l^y — 2/i*£ y z 


E (E l ,E 2 ) [Pa] 


A [N/m] 


A = E 2 S/h 


C (K, G) [Pa] 


A [N/m] 


A = GS/h 




- — l-D s ' Vs 


7] s := 77, v s := v 


i [Pa^s' 1 } 






P[Pa l ~ n ] 


S [N^m- 1 ] 


Ps = ^(V2r 2 P 




7s [N l - n m- 1 ] 


ls:= ^ = ^(V2)"- 2 7 


C [Pa^s- 1 ] 


C s [N^ms' 1 ] 




e = e p + e D 


&s = &s,p ~\~ &s,D 


e SyP = Sh ep, e S) D = Sh eo 



Table 4 



Stress-strain vs. force-displacement rules. In the ideal case of a bearing device subjected to bi- 
directional shear, see Figure 2b, stresses and strains are uniform in the device volume. Both £ s 
and £ are postulated to be of the KBC type. £ s , rj s and v s are expressed as functions of global 
quantities like forces and displacements, while £, 77 and v depend on local quantities. 
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Fig. 1. Classification of Bouc- Wen-type models. 
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Fig. 2. (a) Scheme of the system associated with the 2-DoF model (16). m is the system mass, (b) 
An ideal bearing device subjected to bi-directional shear excitation. 
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Fig. 3. Elasto-damaged BW model, uniaxial stress an = %[dev(a)]n- (a) Strain-history en and £22- 
(b) an stress evolution, (c) Stress-strain loops, (d) The damage evolution with the corresponding 
dissipated energy. 
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Fig. 4. Loading-unloading loop of a BW model, uniaxial stress an = ^[dev(a)]n. (a) Total and 
plastic strains, (b) Stress-strain behavior, (c) The evolution of Y d and of the damage, (d) Dissipated 
energies. 
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Fig. 5. Non-proportional loading of 2-DoF PAW and KBC models, (a) Displacement history for the 
x- and y-directions. (b) Force-displacement loops (u x ,Z x ). (c) Force-displacement loops (u y ,Z y ). (d) 
Damage evolution: the curve is the same for both models; see (19). 
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Fig. 6. 2-DoF Bouc-Wen model with energetic linear degradation rules. Case of stiffness degradation 
proportional to e s : (a) loading-unloading loop; (b) dissipated energies. Note that e SjP is negative; 
(c) damage evolution and strength reduction. Case of stiffness degradation proportional to e S:P : (d) 
loading-unloading loop; (e) dissipated energies; (f) damage evolution and strength reduction. 
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Fig. 7. Pivot rule (thin line) vs. the new stiffness degradation rule (38) (thick line), (a) Displacement 
history, (b) Force-displacement loops for the increasing amplitude phase and (c) for the decreasing 
amplitude phase, (d) Dissipated energy increments for the pivot rule, (e) Damage evolutions for the 
increasing amplitude phase and (f) for the decreasing amplitude phase. 
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Fig. 8. Pivot rule (thin line) vs. the stiffness degradation rule (40), based on the dissipated energy 
by plasticity (thick line), (a) Force-displacement loops for the increasing amplitude phase and (b) 
for the decreasing amplitude phase, (c) Total dissipated energy by the pivot rule and by the plastic 
energy rule. The dotted line indicates the total energy dissipated in the case of the stiffness rule (38), 
shown in Figure 7. (d) Damage evolution. 
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9. App 



Bursietal 



i cation to experimental data, (a) Partial strength beam-to-column steel joint JB1-3A 
2002). (b) Experimental Moment-rotation loops: u=93;, — cpf, Z=Mb eam (c) Predicted 



loops with the rule (38). (d) Predicted loops with the pivot rule (14)2- (e) Damage evolution predicted 
by the model endowed with the rule (38) (thick line) and by the model with the pivot rule (14)2 
(thin line). 
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